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We determine the power spectrum for inflation models covering all regimes from cold (isentropic) 
to warm (nonisentropic) inflation. We work in the context of the stochastic inflation approach, 
which can nicely describe both types of inflationary regimes concomitantly. A throughout analysis 
is carried out to determine the allowed parameter space for simple single held polynomial chaotic 
inflation models that is consistent with the most recent cosmological data from the nine-year Wilkin- 
. . . son Microwave Anisotropy Probe (WMAP) and in conjunction with other observational cosmological 

CO ' sources. We present the results for both the amplitude of the power spectrum, the spectral index 

and for the tensor to scalar curvature perturbation amplitude ratio. We briefly discuss cases when 
, running is present. Despite single field polynomial- type inflaton potential models be strongly disfa- 

04 ' vored, or even be already ruled out in their simplest versions in the case of cold inflation, this is not 

5_( ' the case for nonisentropic inflation models in general (warm inflation in particular), though higher 

^ , order polynomial potentials (higher than quartic order) tend to become less favorable also in this 

case, presenting a much smaller region of parameter space compatible with the recent observational 
cosmological data. Our flndings also remain valid in face of the recently released Planck results. 
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I ^ Microphysics during cosmological inflation is largely believed to provide the key ingredients responsible for seeding 
O . the large-scale structure formation in the Universe. Fluctuations generated during an early phase of inflation provide 
> ' a source of nearly scale invariant density perturbations. This prediction from inflation is found to be nicely consis- 
^ . tent with the cosmological observations, in particular, with those from the measurements of the Cosmic Microwave 

I \- Background (CMB) radiation 

Primordial density perturbations originating during inflation can come from both microscopic flcld dynamics and 
' metric fluctuations. Those from field fluctuations have been identified and previously been studied as coming from two 
^ \ well separated and distinct sources, according to the quantum fleld dynamics involved during inflation. The details 
. of inflation model building, like the fleld content and the couplings of the inflaton field to other fields or to light field 
' degrees of freedom, can largely affect the dynamics during inflation. For instance, they determine how ultimately the 
CO • vacuum energy density stored in the infiaton field end up converted into radiation, allowing the Universe to make the 
transition from the inflationary to the radiation dominated phase. In one extreme case, where the couplings of the 
, inflaton field to other field degrees of freedom can be neglected during inflation (but necessarily become important later 
CO ' on in a preheating/reheating phase), density fluctuations are mostly sourced by quantum fluctuations of the inflaton 
fleld. This is the case of the cold, or isentropic inflation scenario Q. In cold inflation a thermal radiation bath is only 
generated at the end of inflation, after the decay products of the inflaton fleld produced during preheating/reheating 
thermalize. In the opposite case, it may happen that the couplings among the various fields are sufficient strong to 
effectively generate and keep a quasi-equilibrium thermal radiation bath throughout the inflationary phase. In this 
situation, the inflationary phase can smoothly connect to the radiation dominated epoch, without the need, a priori, 
of a separate reheating period. This alternative inflationary scenario was named nonisentropic, or warm inflation (for 
recent reviews, see, e.g., Refs. @|)- In warm inflation the primary source of density fluctuations comes from thermal 
fluctuations, which originate in the radiation bath and are transferred to the inflaton 04^ • 

Thermal fluctuations in warm inflation can be modeled by a stochastic Langcvin equation (SLE), where a dissipative 
term describes the transfer of energy from the inflaton fleld to the thermal radiation bath, while the backreaction 
of the thermal fluctuations in the radiation bath on the inflaton field is described by a stochastic noise term. This 
SLE can be derived from first principles p^| - [l3| and it has been shown p^ - [l6j to appropriately govern the evolution 
of the inflaton perturbations during warm inflation. The formation and maintenance of a thermal radiation bath 
during warm inflation is highly dependent on the microscopic quantum field theory details involved, like the type and 
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magnitude of the interactions of the inflaton field with the other field degrees of freedom and how these other field 
degrees of freedom interact among themselves. 

Specific model realizations of warm inflation and the details of the quantum field theory involved have been reviewed 
recently in 0, Q . Despite the realization of the warm inflation scenario be mostly a model dependent problem, the 
basic conditions for warm inflation dominating over the cold inflation scenario are fairly simple and can be expressed 
in terms of the magnitude of the effective dissipation coefficient T in the effective inflaton's equation of motion, 
the temperature T of the radiation bath (assuming fast enough thermalization of the decay products during warm 
inflation, which is again also a model dependent problem) and on the Hubble scale during inflation, H. Accounting 
for a nonvanishing radiation bath, with energy density pij, an accelerated inflationary epoch in the Universe is only 
possible if pn <C p^, where is the inflaton energy density. Even if is much less than p^, we can still have 

p^^^ > H. But since pa ^ T^, we typically have for warm inflation that T > H. This is the basic condition that 
indicates that thermal fluctuations should dominate over the quantum inflaton field fluctuations, which are Sip ^ H, 
when computing the inflaton power spectrum in the density perturbation calculation. One would think that warm 
inflation would also mean that the dissipation due to radiation production should also dominate over the friction due 
to the metric expansion in the inflaton equation of motion, T > 3H. But this only determines the regime of strong 
dissipation in warm inflation, since in general we can still have dissipation as small as T ~ 10~^H, like for example 
for an inflaton field with a quadratic potential \it\ . or even for values T > 10~^'^H, when including a quartic potential 
for the inflaton and still have a smooth transition from the inflationary phase to the radiation dominated epoch, 
one of the trademarks of warm inflation. 

In nonisentropic inflation models, most notably in warm inflation, there is thus a non negligible contribution from the 
radiation bath to the power spectrum. With the dominant contribution to density perturbations in warm inflation 
models coming from thermal fluctuations. This is opposite to the more common cold inflation models, where the 
quantum fluctuations of the inflaton field make the sole contribution. Little, or no importance, has been given in the 
literature so far to those cases where neither quantum or thermal fluctuations make the dominant contribution to the 
power spectrum, with the possibility of fluctuations from both quantum and thermal sources being equally important 
for the determination of the power spectrum. In those situations, which in fact should cover most of the parameter 
space of realistic quantum field theory models for inflation, we are required to have a complete determination of both 
types of contributions to the power spectrum, coming from quantum and thermal sources, and also to determine under 
which conditions one may eventually dominate over the other. This is a most required study, since, as more precise 
measurements of the CMB radiation anisotropies rapidly become available, greater demand is placed on more precise 
theoretical predictions, such that also better constraints on these models can also be obtained, when contrasting 
them with the observational data. This requires accounting for all important effects contributing to the density 
perturbations. For warm inflation, this requires knowing precisely when the eventual thermalized inflaton fluctuations 
start effectively to dominate over the quantum ones and also to be able to describe all the possible intermediate regimes 
ranging from cold to warm inflation. We here discuss this calculation of the total power spectrum for the inflaton, 
leading to a result that should cover all the regimes ranging from cold to warm inflation. 

As mentioned above, in warm inflation the power spectrum from thermal fluctuations is determined from a SLE 
for the inflaton field fluctuations. For cold inflation there is also an analogou s appro ach. In the stochastic inflation 
program for cold infiation, originally proposed by Starobinsky [l^ (see also |20l-|23j). the quantum inflaton field is 
split in long- wavelength, super-Hubble modes and in short- wavelength, sub-Hubble modes. While the super-Hubble 
component effectively behaves like a classical system, it is influenced by the sub-Hubble quantum modes, constituting 
the environment and that behaves as a stochastic noise term. This leads to a stochastic equation, from which the 
effects of the quantum fluctuations can then be taken into account. This represents an effective coarse-graining 
performed on the quantum inflaton field and on a chosen scale that is at least of the dc Sitter horizon length 1/ H . 

In the warm inflation program, we start by integrating over field degrees of freedom other than the inflaton field. 
The resulting effective equation for the infiaton field turns out to be a Langcvin-likc equation with dissipation and 
stochastic noise terms. We can then further perform now the coarse-graining on the inflaton field according to the 
stochastic infiation program. This way one can describe on equal footing the combined effects of both the infiaton 
quantum and thermal fiuctuations (we will be assuming that a thermalized radiation bath exists, which, as mentioned 
above is a model building problem). By computing the inflaton fleld perturbation power spectrum, we obtain a formula 
that reproduces the cold and warm inflation results as limiting cases, but that can now also describe the intermediate 
regimes between the two scenarios. As an extra bonus of our study, we are able to describe how quantum and thermal 
components of the field fluctuations get entangled (a result that may eventually be relevant to shed light on the 
decoherence problem of cosmological perturbations during inflation, though we will not address this speciflc theme in 
this work). 

This paper is organized as follows. In section |TT] we briefly review the construction leading to the expressions 
describing the dynamics of the inflaton fleld in both the stochastic (cold) inflation and in the nonisentropic (warm) 
inflation cases. In section IIIII the two programs for describing the quantum and thermal fluctuations arc combined 
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and the relevant equations are derived. In section HVl the inflaton field power spectrum for the fluctuations is derived 
and studied. It is briefly discussed the cases when running is present. In section |V] we give our results for simple 
polynomial inflaton field potentials for chaotic inflation (which are the type of potentials that arc presently under 
stronger scrutiny and which are the most constrained by the observational cosmological data) for both the amplitude 
for the power spectrum, the spectral index and for the tensor to scalar curvature perturbation amplitude ratio. Our 
concluding remarks are given in section IVTl An appendix is included where we show some technical details. 



II. LANGEVIN-LIKE EQUATIONS FOR COLD AND WARM INFLATION 



Let us start by briefly reviewing the derivation of the inflaton's equation of motion in the context of the Starobinsky's 
stochastic inflation approach and then, next, for the case of warm inflation. 



A. Stochastic approach for cold inflation 

Assuming that the couplings of the radiation fields to the infiaton field are negligible during inflation, as usually 
assumed in cold infiation, then the equation of motion for the inflaton field is the standard one, 

where <& here is the Heisenberg inflaton field operator, H = a/a and a is the scale factor. We will also restrict ourselves 
to the case of the dynamics during the inflationary regime, in which case the metric is well approximated by the de 
Sitter one, where ds^ = dt^ — a^(f)(dx)^ and a{t) = exp(fft), with the Hubble constant as given by the Friedmann 
equation. 

The stochastic inflation approach is based on a coarse-graining of the inflaton field, where the quantum infiaton 
field is decomposed in a long wavelength part, $>, and in a short wavelength part, $<, 



dV{^) 



= 



(2.1) 



$(x,i)->$>(x,i) + *<(x,t) , (2.2) 

and focus is given on the dynamics of the long-wavelength modes of the field. The short wavelength part of the 
field, $<, summarizes quantum vacuum fiuctuations, representing high momentum modes of the field, with typically 
k ^ kfi ~ aH. This separation of modes is implemented through a suitable filter function. Following the original 
works, e.g. Refs. [l9|-|23|, $< is written in terms of a filter (window) function, which projects out the long wavelength 
modes and it is decomposed as 



$<(x,t) = 4>gix,t) = 



(27r)3/2 



(2.3) 



where (j)]i{t) are the field modes in momentum space, ajj^ and Sk are the creation and annihilation operators, respec- 
tively, and W{k,t) is a suitable window filter function. The filter function W{k,t) in the above equation, in its most 
common implementation, like in the original works on stochastic inflation, is expressed in terms of a sharp momentum 
cutoff, 



W{k,t)^0{k~eaH) , (2.4) 

where e is a suitable small number, such that $<(x, t) refers only to the flelds modes with wavelength much smaller 
than the horizon, k is the comoving coordinate wave- vector and k = |k|. Overall, the results obtained from the 
window filter Eq. (|2.4p . like for the scalar perturbation amplitude and other quantities, are found to be very weakly 
dependent on the value of this parameter e, as we will also explicitly sec in our calculations to be presented in the 
next section. It is also known that a sharp window function, like Eq. (12. 4 p. produces correlation functions that tend 
to decay too slowly at large distances, compared to non-smoothed fields [2J]. The use of a smooth window function 
adds a nonsharp smoothing scale and can give the expected asymptotic behavior for the correlations. For example, a 
Gaussian function (2J] , given by 
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W{k,t) = 1 - cxp 



(2.5) 



has also been used in the hterature. It works nicely as a smooth window function, from which we can recover the 
sharp filter case as an approximation, for modes with wavelengths sufficiently shorter than the horizon (or in the 
approximation e ^ 1). 

In Eq. (|2.3p . the field modes (/)k(i) satisfy 



and 



■ o2 ^1 

— +iH- + e/a^ + {V^^{'^y)) 0k(t)=O, 



<^k(t)<^j;(O-0kW</'k(i) 



(2.6) 



(2.7) 



where {V^cf,cf,{^y)) is an average value of T^00($>) that is taken with respect to the dc Sitter invariant vacuum. In its 
most general form, (j)-\^{t) is given by [25| 



c(0 



aH 



3/2 



^ \ aH 



where h\^\x) is the Hankcl function of the first kind, expressed in terms of the Besscl functions of the first and 
second kind as 



i7«(a;) = J,,(x)+*F^(x) 



(2.9) 



and /.t = ■\/9/4 — m? /H"^^ with -m? = (V^^^ (<&>)) ~ constant. For the high frequency quantum part of the infiaton 
field, and when approximated as a free massless field, = 0, the solution for the modes 0k (0 that follows from 
Eq. (|2.8p . becomes 



H 



(2.10) 



where t = —l/[a{t)H] is the conformal dc Sitter time (ranging from —oo to 0). 

Finally, using the field mode decomposition Eq. (|2.2p in Eq. (j2.1[) . the equation of motion for the long wavelength 
part of the quantum field, $>, becomes 



9$. 



1 



df^ ' dt a- 

where, by keeping only the linear terms in $<, the expression for the term in Eq. (|2.1ip is given by 



(2.11) 



dt a? 

Using the definition Eq. (|2.3p . it can be shown that the comutator for ^q{x,t) is 



(2.12) 



[c,(x,i),e,(x',t')] -0 



(2.13) 



which suggests that S,q behaves like a c-number. It can then be shown that ^q can be interpreted as a (Gaussian white) 
stochastic noise term, where {^q) = and with two-point correlation function (in the vacuum), by using Eqs. 
and (|2.12p . given by 
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(e,(x,t)e,(x',t')> = / -0^e^''-^---'^Rc[Ut)mt')] , (2.14) 
where /k(t) is defined in terms of the filter function W and the field modes in momentum space as 



Equation (|2.1ip can then be interpreted as a stochastic equation for the low-frequency, high wavelength field modes, 
whose effect of the quantum (high-frequency) modes is felt through a stochastic noise term ^g. 

B. Langevin equation for warm inflation 

The derivation of the effective equation of motion for the inflaton field in the warm inflation case also involves a 
selection of a relevant field mode, in which we are interested in the dynamics, while the remaining degrees of freedom 
are taken as environmental ones, similar to the cold stochastic inflation case. The main difference is that now the 
relevant object is the inflaton fleld itself, either the zero mode homogeneous mode, in case we are interested only in the 
background dynamics, or the homogeneous mode plus the purely classical thermal fluctuations, in case of the analysis 
involving the density perturbation from warm infiation. The degrees of freedom that are regarded as environment are 
any other fields in the original model Lagrangian. In short (for a more throughout discussion and extended formal 
details involving the quantum field theory derivation, see, e.g., the reviews [H, Q, or the original papers [26| about the 
warm infiation scenario) , one starts with a model Lagrangian density of the generic form. 



L[<i>,X,Y]=L[^]+L[X]+L[Y]+Lint[<^>,X]+Lint[X,Y] , (2.16) 

where $ is the inflaton field, X are any field or degrees of freedom coupled directly to the inflaton field, while Y can 
be any other fields not necessarily coupled to the inflaton, but are coupled to X. The £int['i', -'f] and £int[-'i^, give 
the interaction among these flelds. The evolution of the inflaton field can be properly determined in the context of 
the in- in, or the Schwinger closed-time path functional formalism j27| . 

In the closed-time path formalism the time integration is along a contour in the complex time plane, going from 
t = — oo to -l-oo (forward branch) and then back to i = — oo (backwards branch). Fields are then identified on each 
of the time branches like, e.g., $i and $2, respectively, and so on for the other fields. Due to the duplication of field 
variables in this formalism, one can now have four two-point Green function that can be constructed with each of these 
fields, with the (l,l)-component identified with the causal Feynman Green function, e.g., Gi^i{x, x') = (T$i(x)$i(a;')), 
where T is the time ordering operator. By functionally integrating over the X and Y fields in Eq. (|2.16p in this 
formalism, we obtain the corresponding effective action for the inflaton field 3>, from which its effective equation of 
motion can be derived Q . This procedure results typically in a nonlocal effective equation of motion for the inflaton 
field <&. This equation displays both dissipation and non-Markovian stochastic noise terms and it is, thus, a generalized 
Langevin- like equation of motion jisj . This is the expected behavior for the produced dynamics of the system (the 
inflaton) , when integrating over the environmental degrees of freedom, here corresponding to the X and Y fields (see 
for instance [27| for a throughout review and discussion about these type of equations and their derivation in the 
context of quantum field theory) . 

It is useful to consider a specific case as an example of this procedure for showing how exactly a dissipation and 
noise terms can emerge in a quantum field problem. More specifically, let us consider a model with a Lagrangian 
density of the form of Eq. (|2.16p . 



+ \ {d,xf \m\X^ - ^$2X2 

- ^i[X.Y] , (2.17) 

where X is a scalar field and L [X, Y] denotes the part of the Lagrangian density containing other fields Y and their 
interaction with X. The Y fields can for example represent radiation fields, whose mass is much smaller than the 
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temperature, my <C T. With Y taken as another scalar field, we can choose a simple form for F], with a 

Yukawa coupling to X, for example, 



Li[X, Y] = hmxXY^ + L [Y] 



(2.18) 



whose specific form for the term involving only Y, will not be relevant for our derivation below. The effective 

action for the $ field is determined by integrating out the fields X and Y. The functional integration over the Y 

(X) 

field will affect the X field two-point Green functions G]j {x,x'), i,j = 1,2, dressing it by corrections coming from 
the Y field. Restricting to corrections to the $ field up to quadratic order in X, the functional integration over 
the X field in Eq. (|2.17p can then be exactly performed. The procedure to obtain the effective action for $ and its 
corresponding effective equation of motion follows analogous steps already shown in many previous references (see, 
e.g., Refs. I, [TMl 

for a more detailed account of this). For the model Lagrangian Eq. (|2.17p . keeping terms up 
to order 0(5^), which is the order where nontrivial results first appear (e.g., dissipation and fluctuation terms), the 
expression for the effective action for $, in the context of the Schwinger closed-time path formalism, becomes: 



r[$i,$2] 



i(a,,<i>i)2-m) 



d'^x 



El 

4 



;9_ 

'16 



(0) [<^l{x)~<^l{x)] 

^1{x)g['^\x - x')Hl{x') - <fl{x)Gif\x - x')Hl{x') 
<^i{x)G\f{x - x'f^lix') + 'fl{x)Gi^\x - x'f^lix') . 



I— I d'^xd'^x 



(2.19) 



At this point, it is convenient to define new field variables: $c = {^i + ^2)/'^ and $a = $1 ^ 'I'2, in terms of which 
the action (|2.19p can be expressed in the form: 



r[$c,$A] = j d^xd^x' 1^^^i^{x)^f[^c]{x,x')^^{x') + ^^{x)0r[^c]{x,x')'^c{x') 



'-^c{x)^a{x) 



4! 



5^{x~x') + Y.^P{x,x') 



where T>p'[^c]{x,x') = ^c{x)^c{x')T}p\x,x'), with 



Tif{x, x') 



o4 r 



FiX){x-x' f --^p'^''\x-x' f 



and 



(2.20) 



(2.21) 



0«[$,](x,x') = [d^ + V"{^,) + I](j"iLai] 5\x x') + E^f )[$J(x,x') , 
where y!'P[^c]{x,x') ^ $c(a;)$c(a;')S^^(a;, x'), with 



(2.22) 



E^f )(x,x') = -^F(^)(a; - x')p'^''\x - x') 0(1 - t') 



(2.23) 



where the functions F'^-^\x,x') and p''-^\x^x') in Eqs. ()2.2ip and p.23|) . are defined in terms of the expectation value 
for the commutator and the anticommutator of the X field, respectively, by 



pW(^,^') = imx),x{x')]) , 

fW(x,x') = \{{X{x),X{x')}) 



(2.24) 
(2.25) 
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and in Eq. (|2.22p . S^j'^^^j = —{g'^/2)G[^\o) is a local correction to the <J> field action coming from the integration 
over the X field. 

The first term in Eq. (j2.20p . the quadratic term in is a purely imaginary term. It can be interpreted as follows. 
We can perform a Hubbard-Stratonovich transformation in the functional partition function generating the effective 
action (|2.20p (see, e.g.. Ref. [13), introducing a stochastic field ^(a;) to decouple the quadratic term in $a, 



{x,x')i{x')+i / d^x^{x)<^>A{x) 



(2.26) 



The effective action expressed in terms of (,{x) has now only real terms. Next, the effective equation of motion for $ 
is obtained from the saddle point equation [lO| 



(5r[$„<j>A,e] 



-0, 



(2.27) 



which, by using the effective action written in terms of ^(x), gives a stochastic equation of motion. 



d^x'Oi^[<^,]{x,x')^,{x')=£,{x) 



(2.28) 



where Oji[^c]{x,x') is defined in Eq. (|2.22p . The equation (|2.28p can be seen as a nonlocal Langevin-like equation of 
motion, where ^(a;) can be interpreted as a Gaussian stochastic noise with the general properties of having zero mean, 
{^{x)) ~ 0, and two-point correlation 



(2.29) 



The equation of motion (j2.28p can be further put in a more suitable form, similar to a non Markovian Langevin 



equation, if we define t!"^\x,x') = Y.'}f^\x,x')9{t — t') and a dissipation kernel T)''-^\x,x') as [iJ 



S](^)(x,x')sgn(t - t') = -^23(^)(a;, x') 



(2.30) 



Then, by absorbing local terms in a renormalized effective potential for $, K;ff.r(<I'), Eq. (j2.28p can now be expressed 
as 



dH, + Kff,r(*c) + j dVDW(x,a;')$c(x') = i{x) . (2.31) 

The dissipation kernel D'^-^\x,x') in the above equation and the noise kernel given by Eq. (j2.26p are related to each 
order, under a space-time Fourier transform, by the fluctuation-dissipation relation 



(p,cj) = 2uj 



n{uj) + - 



(2.32) 



where n{uj) is the Bose-Einstein distribution. 

The resulting nonlocal, non-Markovian equation of motion for Eq. (|2.3ip . can be shown to be well approximated 
by a local Langevin equation of motion when there is a clear separation of timescales in the system, which leads, in 
this local approximation, to a local dissipation coefficient defined by 



(2.33) 
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where, for this specific example we have worked out here, with Lagrangian density (|2.17|) . we have that 



llR[^,]{x,x')=-ig^^l 9{t-t'){[X^x),X^{x')]) . (2.34) 

Complete expressions for T can be found in Ref. [l^ for different interactions and regimes of parameters. 

The validity of the localization approximation for the nonlocal terms in Eq. (|2.3ip and region of parameters where 
this can be achieved in specific models, were discussed at length in many previous references about the warm inflation 
implementation and we will not enter or need these formal details here (for the interested reader, please see, e.g., 
Refs. @, [il] and references there in). In short, if the self-energies terms coming from the X and Y fields introduce a 
response timescale l/F and if $ is slowly varying on the response timescale l/F, 



- « r, (2.35) 

which is typically referred to as the adiabatic approximation, then a simple Taylor expansion of the nonlocal terms 
can be performed. The typical microscopic mechanism for generating radiation in warm inflation is through the 
transfer of energy from $ through X that then can decay into the radiation fields Y . Then, if F > iJ, this radiation 
is thermalized sufficiently fast and, by further requiring that F 3> T/T, it can be maintained in a close to thermal 
equilibrium state. Under these conditions, we can write for the inflaton an effective SLE of motion (dropping the 
subscri pt " c" in the field in Eq. p.3ip and working directly in a Friedmann- Robertson- Walker (FRW) background 
metric [IJl) of the form 



d 



d 



+ (3H + T)---V 



at2 



dt 
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(2.36) 



where Kff,r(^') is a renormalized effective potential for the infiaton (after integrating over the X and Y fields), T is 
the dissipation coefficient and describes thermal (Gaussian and white) noise fiuctuations in the local approximation 
for the effective equation of motion and it is connected to the dissipation coefficient through a Markovian fiuctuation- 
dissipation relation, which in expanding space-time, is given by 



(Ct(x, t)eT(x', t')) = 2TTa^^S{x - x')S{t - t') , (2.37) 

where the average here is to be interpreted as been taken over the statistical ensemble. 

The dissipation coefficient T in Eq. (j2.36p depends on the specifics of the model Lagrangian parameters entering in 
Eq. (|2.16|) and can be a function of temperature and inflaton field. Recently, it was derived in Ref. [l^ the different 
dissipation coefficients arising from many different types of couplings and for different temperature regimes that can 
be found in warm inflation, with other relevant regimes for warm inflation also identified in Ref. [29j . For example, for 
the models treated in the Refs. 0,[23, the dissipation coefficient T was found to have the following generic dependence 
on the inflaton field amplitude 0, temperature T and on the mass mx for the X fields coupled to the inflaton, 

= C</'^?7r. c + 2a-26=l, (2.38) 

where the values of the constants a, 6, c depend on the specifics of the model construction for warm infiation and 
on the temperature regime of the thermal bath. For example, for a trilinear interaction of the inflaton field with a 
bosonic X field and at low temperature regimes (for the temperature taken with respect to the intermediate field X in 
the Lagrangian density model Eq. p.l6|) . T < mx), it is found that a = c = 3, 6 = 4, while for a quadratic coupling 
of the inflaton field with a bosonic X field like the one taken in the example shown above and at high temperatures 
(T > mx), it is found that a = l,& = 0, c = — 1. Other types of dependence of the dissipation coefficient on the 
temperature and on the inflaton amplitude can be found for other forms of couplings and regimes of parameters. 

From a model building perspective, warm inflation has several advantages. For example, we can implement La- 
grangian densities of the form (|2.16p in terms of supersymmetry (SUSY) 0, 0, H^l , so couplings of the infiaton field 
with the other (bosonic and fermionic) fields X and between X and Y can be made fairly large, 0(0.1), and we can 
still keep quantum corrections to the infiaton potential small using SUSY. Furthermore, temperature corrections to 
the inflaton potential can as well be made small when T < mx- 
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III. STOCHASTIC APPROACH FOR COLD AND WARM INFLATION 



We now want to combine both the Starobinsky stochastic inflation approach and the warm inflation scenario 
in a single SLE description. Note that this could be implemented from the very beginning in the warm inflation 
construction of Eq. (|2.36p . where besides of integrating over the X and Y fields, a functional integration of the 
sub-horizon modes of the inflaton field $ could also be performed, leaving only an effective equation describing the 
super- horizon part of <&. A functional integration of the sub- horizon modes of the inflaton (in the absence of any other 
fields coupled to the infiaton) and the construction of the influence functional for the super-horizon modes, in the 
context of the cold inflation only, has been done, for example, in [sTI - fssj . This functional integration procedure, when 
applied to cold inflation, is able to describe the inflaton dynamics in a more rigorous way than the simple field split 
implemented in the Starobinsky program and it has been shown to lead to additional contributions in the infiaton 
equation of motion for the long wavelength modes, Eq. (|2.11|) . However, it has also been shown that these extra 
contributions only produce small effects and they can, therefore, be neglected [3l|, [3^. Therefore, the result given by 
Eq. (|2.1ip can still be used. 

By describing the influence of the thermal radiation bath on the inflaton dynamics, according to the SLE (|2.36|) . we 
can now account concomitantly also for the effect of the sub-horizon quantum fluctuations on the dynamics for the 
long wavelength modes of the inflaton field. Hence, the infiaton field $ in Eq. (|2.36p is again split like in Eq. (|2.2p . 
in a sub-horizon, quantum fluctuation part, (/)g(x, i), which can again be described by Eq. (|2.3p . and in a super- 
horizon part, that behaves essentially as the classical inflaton mode. The classical inflaton mode is in turn split in 
the homogeneous inflaton field (f>{t) and in the fluctuations S(fi{x,t). The background homogeneous inflaton can 
be seen as the coarse-grained infiaton field over approximately all de Sitter horizon size n = l/H, 



m = ^ lj'^H^,t) , (3.1) 

and S(p{x,t) are the classical fluctuations on top of it. 

In the absence of the thermal backreaction from the radiation bath, or when T <g; iJ, it is the quantum noise 
that is expected to predominantly contribute to Stpix^ t). Here, we will not attempt to study how exactly the quantum 
fiuctuations turn into classical ones, the decoherence problem of the fluctuations in cold inflation, but assume that the 
decoherence time is short enough that when the inflaton fluctuations cross the horizon they behave or are essentially 
classical fluctuations already, as it is generically expected (the decoherence problem in inflation was also extensively 
studied in the context of stochastic inflation and the influence functional method. For some references, see, e.g., 
prl [31I [3^ and references there in). In addition, as the temperature of the radiation bath increases, or when the 
coupling of the inflaton flcld with the radiation bath becomes sufficiently strong, we expect that now the backreaction 
of the thermal radiation, represented by the thermal stochastic fiuctuations in Eq. (j2.36p . will become more and 
more important, till it eventually surmounts that from the quantum stochastic fiuctuations. Next, we will study how 
exactly these two effects compete in giving the dominant contribution to the infiaton fluctuations (5iy9(x, t). Towards 
this goal, and as explained above, we can then write the inflaton fleld as 



$(x, t) = (j){t) + Sifiix, t) + (f)g{x, t) . (3.2) 

By substituting Eq. (j3.2p in Eq. (|2.36p . and expanding in both perturbations 5ip{x,t) and c/jqlx^t) to first order, we 
obtain the equations for the background field (j){t) and for the fiuctuations dip{x.,t), given, respectively, by 



^ + [3H + Tm^ + V4^)^0, (3.3) 

^ + [3H + Ti^)] ^ ~ + TA^)^ + V^'j Sip = iq + , (3.4) 

where is the quantum noise term Eq. p.l2p , but including now also the contributions coming from the dissipation 
coefficient due to the integration over any of the radiation bath fields (for example, the fields X and Y in Eq. (|2.36p ). 

^ + m + T(0)] ^ - + T,40)0 + V^^iA <Pq . (3.5) 



We note that we will be working always in the local approximation for the SLE, in which case, both noise terms 
in Eq. (|3.4p arc Markovian stochastic noises. We should also note that in this case, a quantum noise term that we 
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would expect to be produced alongside the thermal noise ^t, from the functional integration of the the fields X and 
Y , vanishes identically in this local approximation [isj and the quantum noise term comes exclusively from the 
inflaton field quantum modes, analogously as in the standard stochastic inflation program (for a study of the effects 
on the inflation dynamics of a quantum noise in the non-Markovian case see, e.g., Refs. (ssl. [36j). 

In terms of the modified quantum noise ^q, it can be verified that the commutator (|2.13p still satisfies 



4(x,i),4(x',t') =0 



(3.6) 



which again is suggestive of interpreting as a classical quantity. 

From the equation of motion for the background field 0(t), Eq. p.3p . we have the slow- roll equation as that of 
warm infiation, 



iH{l + Q)^ + V^{(j))=Q , 
where we have defined Q = T/3H, and we can also define the slow-roll coefficients 



(3.7) 



1 



16nG 



Ya 

V 



< 1 + Q 



(3.8) 



SttG V 



< 1 + Q 



and 



(3.9) 



(3.10) 



where the third slow-roll coeflacient, /3, is typical of warm infiation, accounting for the possible dependence of the 
dissipation coefficient on the inflaton's amplitude. 

Going to momentum space, it becomes convenient to express the equation for the fluctuations, Eq. p.4[) . in term 
of the variable z = k/ (aH). If we also use the definitions of the slow- roll coefficients, Eqs. p.8p . p.9p and p.lOp . the 
equation satisfied by the fluctuations Sf(k., z), the Eq. p.4p in momentum space and in terms of the variable z, can 
be written in the form 



V'(k,z)--(3Q + 2)V(k,z) + 



eT(k,z)+e,(k,z) 



1+3 



r,-pQ/il + Q) 



Sipik, z) 



(3.11) 



where primes in </3(k, z) mean now derivatives with respect to the variable z. Equation p. lip for the inflaton fluctu- 
ations generalizes the one obtained in [isj . Its general solution can be expressed in terms of a Green function. 



(5(p(k,z) = / dz'G{z,z') 



1^2 



i^iz')+^T{z') 



where the Green function G{z, z') is given in terms of Bessel functions. 



G(z, z') = -z-'z"' [Mz)Y^{z') - Mz')Y^{z)] , 



(3.12) 



(3.13) 



with z' > z and 



i/ = 3(l + Q)/2, 

/ , ^ 3/3g 

a = \ 1^- + — 

V i + Q 



(3.14) 
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Using Eq. (|3.12p . we can write the following expression for the correlation {SipCk' , z)5ip{'k' , z)), 

{6ip{k, z)6ip{k' , z)) = — / dz2 I dziG{z, zi)G{z, Z2) 5 Z2 (C« (k, ^i)^g(k', 2:2)) 



1 
IF 



dzo 



dziG{z, zi)G{z^ Z2 



^1 ^2 



-(eT(k,ziKT(k',z2)) . (3.15) 



Note that there is no inter-correlation between the two stochastic noise terms, the thermal one S^t ^md the quantum 
one ^g, since they come from different sources, as we have already discussed before. Below, we evaluate the two terms 
in the right- hand- side of Eq. p.isp separately. 

IV. TOTAL POWER SPECTRUM INCLUDING QUANTUM AND THERMAL NOISES 

Let us now compute in details each one of the contributions to the power spectrum appearing in Eq. (j3.15p . 

A. The quantum noise component for the power spectrum 



Let us now evaluate the quantum noise component for the power spectrum, which comes from the first term in 
right- hand-side of Eq. (|3.15p . Here, we will follow a similar approach as used in [s^l (see also [l^l) for the computation 
of the power spectrum in the stochastic cold inflation case. We first need to compute the quantum noise two-point 
correlation, (^^(k, zi)^q(k', Z2)). This is determined by the equation defining the quantum noise ^g, Eq. p.Sp . 

From the Fourier transform of Eq. p.Sp and working in terms of the comoving time r = — 1/ {aH) , which significantly 
facilitates the calculations, we find 



4(k,r) = -i/V 



Q^2 ^ 5^ + + ^2^2 + H2^2 



,(k,r) 



(4.1) 



where 



(^^(k, t) = W{k, t) Uk(T)a„k + <^k(T)4 



(4.2) 



Under the assum pti on that the quantum fluctuations are [22j 0{-\/h) and also considering that the dissipation term is 
at least of order [1J| 0(?i), then the quantum field modes t/ik still satisfy the usual differential equation (in conformal 
time) at leading order: 



a2 



2 d 



with solution of the form of Eq. (|2.8p . which in conformal time is 



Mr) = , 



(4.3) 



Mr) = ^{\r\r/'Hi'Hk\r\) , 



(4.4) 



where ^ = \/9/4- 3?? « 3/2 - -q. 

From Eqs. ()4.ip and ()4.2p . we obtain for the two-point correlation for the quantum noise the expression 



{i,{k,r)Uk',T')) = {TT'rH' f/k(r)/^,(T')(a-k4') + /k(r)/k'(r')(4'a-k) 



(4.5) 



where fk{T) is defined as 
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Mr) 



W" 



2 + 3Q 



w 



3/?Q 



(1 + Q)t- 



,w 



0k(r)+2W^yk(r) . 



(4.6) 



In Eq. (j4.5p . the averages involving the creation and destruction operators, (a-kOk/) and (aj^,a_k), are evaluated in 
a thermal radiation bath, so they are in fact quantum statistical averages. We use the definitions, (a-kOk/) = [n{k) + 
l](27r)^(5(k + k') and (aj^,a_k) = n(fc)(27r)'^(5(k + k'), where n{k) is the distribution function for the high-frequency 
quantum inflaton modes. In warm inflation, it assumes the thermal Bose distribution form, n{k) = l/[exp(/c/(aT)) — 1]. 
It is been implicitly assumed here that the inflaton fluctuation momentum modes crossing the horizon have had enough 
time to reach thermal equilibrium with the radiation bath. This condition of thermal equilibrium for the inflaton 
fluctuations is mostly a model dependent problem, requiring the study of the dynamical microphysics involved, the 
details of the interactions involved and the knowledge of the held content and properties of the radiation bath. 
However, it can be verifled and checked that thermal equilibrium can be achieved in general through the study of 
the Boltzmann dynamics for the particle's distributions, which can be obtained from the determination of the rate of 
energy transfer from the background inflaton field to the radiation bath and whose formal details can be found, e.g., 
in Ref. [11] and also discussed in Ref. [29| . 

The expression (j4.5|) for the two-point correlation for the quantum noise can then be written as 



(4(k, r)C;(k', r')) = (2^)35(k + ]^){TT'fH' [2n{k) + 1] Re [/k(r)/^, (r') 



(4.7) 



which can also be expressed in terms of the variable z = k/{aH). Using the definition of the inflaton power spectrum 
Ps.^ in terms of the two-point correlation function for the fluctuations as Q 



k^ f (Pk' 
P^v = ^J (^('5<P(k,^)'5<p(k',^)) , 

we then obtain for the quantum noise contribution for the power spectrum the result: 



(4.8) 



6ip 



fc3 



d^k' 



2Tr^H^ J {2tt) 
[2nik) + l]^\F,i 



dz2 



dziG{z, zi)G{z, Z2){ziy 



'i-2^(z2)-^-'''(e;(k,zi)e;(k',z2)) 



(4.9) 



where 



Fkiz) ^ J dz'G{z,z'){z'f-^'' 
+ 2I^,,,(z')0k,,'(^')} , 



(4.10) 



and G(z, z') is given by Eq. (f3T3| . 

By integrating by parts the term with W^z'z' in Eq- (|4.10p using the identity for the derivative of the Bessel functions 

(where Za ^ Ja,Ya,H2^^), 



= Za-l{z) Za[z) 

dz z 



(4.11) 



and using Eq. ()4.4|) . after some algebra, we obtain for Eq. (|4.10p the complete result in terms of an arbitrary filter 
fimction W{z), 
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-3/2 VTJ roc 

^'^W = 4fc3/2 J dz'iz'f/'-- [Mz)Y^{z') - Uz')Y^{z)] 
J9^W,A^')Hl'\z') + z'W,Az')Hl'Uz') 



3Q ( pQ 



z' \l + Q 
n^/^z'^H 



+ v] W{z')Hl'\z') - 3gW^(z')<2i(^') 



4fc3/2 



dz'{z'f'^-Vdz)Yo.-i{z') - J^^,{z')Y^{z)] W^,,{z')H^^\z'). (4.12) 



If wc now make use of the window function (j2.4p . the integrals in Eq. (|4.12p can all be performed. If we also neglect 
terms proportional to the slow-roll coefficients and taking the leading order terms in the resulting expression, using 
that 0<z<e<l, we obtain for Eq. ()4.12p the much simpler result 

F.(.)--^^[l + 0(.^)] . (4.13) 

We note the important result that at leading order Fk{z) is independent on the dissipation coefficient Q. Dependence 
on Q in Eq. (|4.13p can be verified that only appears at order 0(6^) onwards and it involves, multiplicatively, the 
slow-roll coefficients, being, therefore, negligible those contributions. Consequently, the quantum component of the 
spectrum results to be only very weakly dependent on the dissipation. 

Finally, from Eqs. (j4.9p and (j4.13p . the quantum noise contribution for the power spectrum then becomes 

which, as conventional, it is evaluated at the moment of horizon crossing, z = k/{aH) = 1. Note that the thermal 
enhancement factor in the above result was also found by the authors in [39| . although in there it was assumed 
already a thermal initial distribution before inflation and that the thermal background was assumed decoupled from 
the inflaton dynamics (i.e., no nonisentropic or warm inflation regime). Here, we are only requiring a thermal radiation 
bath during inflation (as is generically the case in warm inflation models) and that the inflaton fluctuations are in 
thermal equilibrium with the radiation bath at the time where the relevant scales leave the horizon, as already 
discussed before (for a discussion of the different microscopic processes possibly contributing and leading to this 
thermal state, see, e.g., Ref. [29l|). 

Using that 3 — 2^ ~ 2?7 ^ 1 and taking T = in Eq. (|4.14p . we obtain the standard, approximately scale- invariant, 
spectrum result for cold inflation [3j 

P^T' ■ (4-15) 

If we choose another window function, for example the Gaussian one Eq. (j2.5p . we have to perform the integrals 
that appear in Eq. (|4.12p numerically. We have explicitly checked numerically that for values of the parameter e < 0.1 
in the Gaussian filter Eq. (|2.5p . the results obtained for the power spectrum are in very good agreement with those 
obtained with the step function filter, Eq. ()4.14p . In fact, when e < 0.1, the step filter is a very good approximation to 
the Gaussian filter. For example, this can be shown from the numerical results for the normalized spectrum obtained 
with the Gaussian filter, with the thermal contribution factor out, and whose results are shown in the figure [TJ 



B. The thermal noise component for the power spectrum 

Let us now compute the contribution in Eq. (j3.15p depending on the thermal noise correlation function, given by 
the second term on the right-hand-side of Eq. p.lSp . First note that the thermal noise correlation function Eq. ()2.37p . 
in momentum space and in terms of the variable z, can be written as 



(|T(k, 2i)|T(k', Z2)) = 2TT^zlz25 {z, - z^) {2nfS{k + k') , 



(4.16) 
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FIG. 1: Normalized power spectrum obtained using a Gaussian window filter, with e — 0.5 (solid line), 0.3 (dashed line) and 0.1 
(dotted line). 



where we have used the time variable written in terms of z, 



l/H 



(4.17) 



together with 



S{ti - t2) = H5 
and the property of the Dirac delta-function, 



1 f 



(4.18) 



where xq are the roots of f{x). 

From Eq. (|4TT6)) and using the definition of the power speetrum Eq. (|4.8p . we then obtain the contribution to the 
power spectrum coming from the thermal noise source as given by 



^Sv = 27r2g4 y (^:;;^ y ^^^2 y dziG(z,zi)G(z,Z2) ^ ^ (^T(k,zi)CT(k ,22)) 

^ ±£ / dz'z'2-4-G(z,z')' . (4.20) 

Moss and Graham in Ref. [l^ have also found a similar result for the power spectrum, except that they have, 
instead of T in their result, an effective dissipation coefhcient Toff. They then determine this effective dissipation 
coefficient by a matching procedure with the expected spectrum obtained for a flat space thermal field theory. Here, 
as we are going to see below, this is not necessary, since we are already determining the complete inflaton power 
spectrum, that includes both thermal and quantum noise sources. 

For the relevant observational scale for the inflaton perturbations, which are produced in the late stages of inflation, 
the value of the argument z in (|4.20p can be taken as been small. Then, for example, we can approximate the integral 
appearing in Eq. (|4.20p such that 
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/•CXI 2 pca 

Jz 4 Jq 



'{z'f-^'^Jliz') 

= - [^^y.{^^)? ^^r'^''^W^"'\. , (4.21) 

4^ '^^ 20Fr(i^- l/2)r(a + i^- 1/2) ' ^ ' 

where T(x) is the Gamma function. Likewise, we can also use the asymptotic form for the Bessel function, i.e., 

By neglecting terms proportional to the slow- roll coefficients, with 77, /3 <§; 1, such that we can take a ~ = 
3(1 + Q)/2 in the arguments of the Gamma functions in Eqs. (|4.2ip and (|4.22p . the thermal component for the power 
spectrum, Eq. (|4.20p . can finally be approximated as 



~ 16^"^ r(2i^- i/2)r(j.- 1/2) 



2r,-2/3Q/(l+Q) 



HT 3Q8'^ [r(3/2 + 3g/2) 



l3 



(4.23) 



7r2 (l + 3Q)r(5/2 + 3Q)r(l + 3g/2)' 
Note that the thermal spectrum is still approximately scale- invariant, as we would expect. 

C. The total power spectrum and relevant cosmological parameter results from observations 

From the result for the quantum contribution to the power spectrum Eq. (|4.14p and the result for the thermal 
contribution, Eq. (|4.2ip with Eq. (|4.22p . our result for the total power spectrum becomes, 



, , HT 



3Q_.2a^,2.-2c r {af T{v-l)T{a-v + 3/2) ^ H ^^^^ f zH , 



T{i^-^)r {a + 1^-1/2) T \2T 



(4.24) 



which is conventionally evaluated at the moment of horizon crossing, z = k/{aH) = 1. Note that in the zero 
temperature limit, T — > 0, we also have that [l^ T (Q — ?> 0) and the cold inflation result, Eq. (|4.15p is re- 
obtained as expected. In the limit of Q ^ 1 and T 3> -ff , we obtain Psip oc HT, which is the original result obtained 



for warm inflation by Berera and Fang in Ref. 126L In the limit Q ^ 1 and T ^ 1 we obtain that Pg^p oc TViJT, 
which is the result found by Berera and Taylor in 3] and also by Hall, Moss and Berera in Q . 

From a perturbed Friedmann- Robertson- Walker metric with spatial curvature perturbation -tp, it can be defined 
a gauge invariant comoving curvature perturbation as !R(x, t) = 'ip{x,t) + H5ip{x.,t)/(j), which remains conserved 
throughout the evolution on super horizon scales. We can then choose a gauge, for example a spatially flat gauge, 
and write the curvature perturbation simple as Q 3^ = HSip/4>- Thus, from Eq. (|4.24p . we obtain the amplitude for 
the curvature perturbation during inflation. Ax, as 

- -^Ps^ = A^(fco) (^^j , (4.25) 
where ko is a pivot scale [1], fco = 0.002 Mpc~^, is the spectral index, defined by 

dlnA|_dlnA| 

where Nf, denotes the number of e-folds of inflation. These are all quantities that are read at the longest length 
scales of interest and taken to be those that first crossed the Hubble radius and which corresponds to the scales that 
crossed the horizon iV* e-folds before the end of inflation. N^, is between the usual canonical values, N^, = 50 and 60 
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e-folds (40[ , whose specific value depends on the details of inflation and on the reheating phase after inflation. These 
will be the values of e-folds before the end of inflation that we will consider in the numerical analysis to be shown in 
section M below. 

The most recent CMB radiation measurements from the WMAP nine years of accumulated data (WMAP9yr) 
gives for A|,(fco) and n^, respectively, the values Q: A|(A:o) = (2.41 ± 0.10) x 10^^ and = 0.972 ± 0.013, 
for the measurements from WMAP only. When the results from WMAP are combined with the results obtained 
from the extended CMB data (cCMB) (the CMB data coming from the Atacama Cosmology Telescope (ACT) and 
the South Pole Telescope (SPT)), from baryon acoustic oscillation (BAO) results and measurements for Hq from 
galaxies surveys, it is obtained (assuming also contributions from tensor modes) A^(fco) = (2.464 ± 0.072) x 10~^ 
and Ug = 0.9636 ± 0.0084, respectively, which slight improves over the previous WMAP seven years of accumulated 
results (WMAP7yr), combined with similar databases (WMAP-|-BAO-|-i7o)i which have produced the results Q 
A|(fco) = (2.430 ± 0.091) x lO'^ and = 0.968 ± 0.012. 

Another cosmological parameter that helps to constrain many inflation type of models is the tensor to scalar 
curvature perturbation ratio, r. 



4^ , (4-27) 



- A^(fco) 

where A^(fc) is the power spectrum of tensor metric perturbations, which for cold inflation is given by Q 



where Mp = l/y/SnG is the reduced Planck mass. The value of r is still poorly known, but an upper limit can be 
found from the observational data. The recent WMAP9yr results (WMAP+eCMB+BAO+Fo) gives r < 0.13 (at 
95% CL). 

It has been shown in Ref. [4l| that in a thermal environment there can be a contribution to the tensor perturba- 
tion (|4.28p coming from stimulated emission, if gravitons were in thermal equilibrium with the radiation bath. This 
effect would then imply a thermal enhancement factor coth [ziJ/(2T)] in Eq. (|4.28|) . just like the same factor also 
appearing in the quantum scalar perturbations, Eq. (|4.14p . The presence of this thermal enhancement factor in the 
tensor perturbation power spectrum depends whether the gravitational wave field is in thermal equilibrium with the 
radiation bath or not. In the original assumption made by Bhattacharya, Mohanty and Nautiyal in Ref. jZlj . inflation 
was considered to start in a state that was initially in a radiation dominated epoch, with temperature T. Alternatively, 
if a thermal radiation bath develops in the course of inflation (as assumed in warm inflation) and the gravitons have 
enough time to equilibrate with the radiation bath and be strongly enough coupled to the particles in the radiation 
bath, then we also expect them to develop a thermal radiation-like spectrum, though in this case we would typically 
expect to be necessary at least temperatures for the radiation bath of order [i^ T ~ Mp, which may be too high 
compared to the typical temperature scales during warm inflation, which are of the order of the Grand Unification 
(GUT) scale and below it. Nevertheless, the presence of the thermal enhancement factor coth [z H / {2T)] appearing 
multiplicatively in Eq. (|4.28|) would increase the tensor to scalar curvature perturbation ratio and can potentially 
produce a much stronger constraint on nonisentropic inflation models than in the absence of this factor. Since some 
authors working on warm inflation have included explicitly in their studies this factor in the tensor curvature spec- 
trum (see, for example, Refs. [H, I3) and given the status of the presence of this term is unclear in the context of 
nonisentropic inflation in general, we will consider cases with and without the thermal enhancement factor included in 
the tensor curvature spectrum in the results to be presented below in section |Vl so both situations can be compared 
with the present observational constraint on r. 

The specialized expression for r in the context of stochastic nonisentropic inflation model, with the thermal en- 
hancement factor included in Eq. (j4.28p (and with the spectra taken at the moment our current horizon scale crossed 
the dc Sitter horizon during inflation, i.e., = k^/{a^H^) = 1) is then given by 



(1 + Q^n^Ps, 



■ coth 



(4.29) 



where wc have used the expression for the amplitude of scalar curvature perturbation Eq. (|4.25|) . along with the 
slow- roll equation for (p coming from Eq. p.7p and the equation dcflning the slow-roll coefficient e, Eq. 
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It is useful to check some of the limiting behaviors for r with respect to the dissipation and temperature. The result 
)r r becomes m( 
(at z = = 1), 



for r becomes more transparent if we instead look at the inverse ratio l/r. Using that Pgtp — -Pj^"^ + Pf'^\ we obtain 



1 (1 + Q)^ + 



r 16e 4£iJ2coth(^) 



(4.30) 



Taking for example the weak dissipation regime Q <C 1 and for either T <^ H or T ^ H, wc obtain l/r = 
l/(16e) + 0{Q^), i.e., r ~ 16e, which is the standard cold inflation result For Q ^ 1 and T ^ H, wc obtain that 
l/r = •\/37r(3^/2/(32e) + 0{Q^), i.e., r ~ case of not having the thermal enhancement 

factor in Eq. (|4.29p . or also when T <C -ff in (j4.29p . wc obtain the warm inflation result for r in the strong dissipation 
regime, r ~ 16eiJ/(y3^rQ5/2). We can then conclude that even in the presence of the thermal enhancement factor in 
r, for the weak dissipation regime Q <C 1 the cold inflation result for r is preserved independently of the value of the 
temperature for the thermal radiation bath. It is only in the strong dissipation regime Q 3> 1 that the ratio is affected 
by the dissipation, which makes r smaller by a factor • This leads to an important aspect of warm inflation, 
making it possible to evade the so called Lyth bound [43 (in particular, as observed in Rcf. [iBi, this implies that 
large field inflation models do not necessarily yield a large tensor-to-scalar ratio). Therefore, thermal noises only affect 
scalar type curvature perturbations, but primordial gravitational waves (and associated tensor perturbations) are still 
originated from quantum vacuum fluctuations (though it may receive a correction dependent on the temperature, 
when in a thermal bath, because now the vacuum becomes a quantum statistical vacuum at some temperature T and 
the gravitons can have a non zero occupation number, as shown in Ref. [4l|). 

For the sake of completeness, we also show here the limiting behaviors of rig in relation to Q and T. Details 
concerning the derivation of ti^, in the context of nonisentropic inflation in general, are given in the appendix A. From 
Eq. (|4.26p . taking the Q = and T = limit, we obtain the expected result from cold inflation 

n,s = 1 + 2?7 - 6£ . (4.31) 

In the large Q and T limit (and temperature independent dissipation), we recover the standard warm inflation result 
(e.g., from Hall, Moss and Berera Ref. Q): 



1 



9 9,^ 3 ■ 

£ B + -n 

4 4'^ 2 ' 



0(1/Q^-^) + 0(i7/T) . (4.32) 
Finally, in the low dissipation Q ^ 1 limit, we obtain the first order correction for Ug: 



= 1 + 277 - 6£ 4- (8£ - 2r,)Q + 0(0^) . (4.33) 

If we relax on the assumption of a power spectrum that is a pure power-law, like the one given in Eq. (|4.25p . and 
one admits a running spectral index, then wc can write j47| 



^ ns{ko) — l+ ^ \n{k/kQ)dns /d\n k 



Ai{k) = Ai{ko)i^-j . (4.34) 

With the possibility of a running, the spectral index tends to be blue tilted, as opposite to a red tilt {ug < 1) in the case 
of no running. Some inflation models, like hybrid inflation, in fact prefer a blue tilted spectral index, while others, like 
chaotic inflation type of models, tend to be more red tilted. Observational data, however, arc still not precise enough 
to give a definitive answer whether a running is favorable or not, with no statistically significant deviation from a pure 
power-law spectrum. For example, the WMAP9yr results for the running are only with less than 3a significance 0]. 
Tensor mode amplitude and scalar running index parameter, when fitted together, give Q: r < 0.47 (at 95% CL), 
dus/dlnk ^ -0.040 ± 0.016 and = 1.075 ± 0.046 (using the combined data from WMAP9yr-HeCMB+BAO+i/o). 
For the analysis and results to be presented in the next section, we have chosen to concentrate on the cases of no 
running of the spectrum, when contrasting the results obtained for nonisentropic inflation with the observational data. 
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V. RESULTS AND CONSTRAINTS FOR NONISENTROPIC INFLATION MODELS FROM THE 

OBSERVATIONAL DATA 

Since the early releases of the WMAP CMBR results, a significant pressure has been put on inflation models. Both 
the spectral index and the tensor to scalar perturbation ratio results have shown to be able to exclude large classes 
of inflationary models. For example, hybrid inflation models are disfavored for predicting in general a too blue tilted 
spectrum, while simple chaotic polynomial type of potentials, in the case of cold inflation, even though they have a 
preference for a red tilted spectrum, they have been seen to poorly fit the observational data (though nonminimal 
coupling to gravity can make them compatible with the observational data Q). In particular, the recent WMAP9yr 
results already exclude (at the 68% CL) all simple chaotic cold inflation polynomial potentials V{(f>) oc <f>P, with power 
p > 1 (see, for instance, figure 7 in Ref. Q). 

Since most of the present constraints from the observational data have been applied most to cold inflation models, it 
is also important to have a throughout analysis of how general nonisentropic models withstand the recent observational 
data as well. Even though some previous works (see, e.g., Refs. [1, |H, 13) have addressed the consistency of warm 
inflation with the observational data, only a few examples with very specific values of parameters were studied. A 
complete analysis covering larger regions of parameter space and ranging from cold to warm inflation is still missing. 
Since the larger pressure from the observational date are on chaotic polynomial type of inflation models, which by 
themselves represent an important class of inflaton potential models from a microscopic model building perspective, 
these arc, therefore, the type of inflation models we will be interested in studying next. 

In the following wc will study nonisentropic chaotic inflation potential models with a generic potential of the form 



where A is a dimcnsionless parameter flxing the energy scale and we will consider the cases p = 6, 4 and 2, which are 
the most common type of chaotic inflation potentials used, for example, in cold inflation studies in general. 

Having fixed the class of infiaton potentials that we will analyze, we next need also to specify the dissipation 
coefficient T appearing in the expressions. Nonisentropic inflation can be regarded as a two-fluid system when 
evaluating density perturbations, where both the inflaton and the radiation fluctuations are effectively coupled through 



the temperature dependence of the dissipation coefficient T. It has then been shown in Ref. [15| that because of 
this dependence of the dissipation coefficient on the temperature, the radiation perturbations can back react on the 
infiaton perturbations, generating undesirable growth modes. However, recently in Ref. [l^, it was found that intrinsic 
microscopic decay processes in the radiation bath itself can cause it to slight depart from equilibrium, giving rise to 
viscous dissipative effects in the radiation fiuid itself. These viscous effects were then shown to be able to remove the 
growing modes and to lead to a power spectrum for the perturbations consistent with cases where the temperature was 
absent in T. In this case, where the temperature can be neglected or it is absent in the expression for the dissipation 
coefficient, there are two representative forms of dissipation coefficient that can emerge and that are consistent with 
microscopic derivations (see, e.g., Refs. [13, 113) • ^ constant dissipation coefficient T = Tq and a dissipation that has 
a quadratic dependence on the infiaton field, T oc (j)"^ (these two forms of dissipation coefficients were also considered 
by the authors of the Refs. [i^, in their studies involving warm inflation). These are, therefore, the two forms of 
dissipation coefficients we will use in our analysis, a constant dissipation coefficient, T = Tq and also one proportional 
to the square of the inflaton amplitude, T = C^(tP' /Mp, with constant coefficient C^. 



A. Cold and warm inflation dominated regions 



As a preliminary study, let us first establish which regions of parameter space, in the plane defined by the dissipation 
coefficient over the Hubble parameter, Q = T/{3H), and by the temperature over the Hubble parameter, T/H, the 
spectrum is dominated by quantum fluctuations or by thermal fluctuations. It has generically been accepted that 
warm inflation is determined by the parameters for which T > H, considered the condition for which the thermal 
fluctuations dominate the spectrum. Warm inflation has also been considered in both the weak dissipation regime, 
(5 <C 1 and in the strong dissipation regime, Q ^ 1. However, it has never been exactly established for which values 
of parameters we could really discriminate a cold inflation like spectrum, where quantum fluctuations dominate over 
the thermal ones, and the opposite regime of warm inflation. 

We have determined in the previous section the contributions to the power spectrum coming from either quantum 
or thermal sources, with relevant results given by Eqs. (|4.14[) and (|4.23p . respectively. We can, thus, now determine 
the region of parameters where one contribution to the spectrum dominates over the other. In figure [5] we show the 
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region of values of Q and T /H for which -Pj^" > , i.e., for which quantum fluctuations dominate over the thermal 

fluctuations, and for the opposite case, P'j^"'' < Pg'^\ i-e., for which thermal fluctuations dominate over the quantum 
fluctuations. 




Q 

FIG. 2: Parameter space, in terms of the temperature ratio T/H and dissipation ratio Q = T/(3-ff), for quantum fluctuation 
dominated spectrum (cold inflation) and for a thermal fluctuation dominated spectrum (warm inflation). Red lines are for 
a constant dissipation coefHcient, while blue lines are for a dissipation coefficient dependent on the square of the inflaton 
amplitude. The dashed lines are for p = 2, solid lines are for p = A and dotted lines are for p = 6 (a quadratic, quartic and 
sextic inflaton potentials, respectively). 

We believe that figure [2] is very illustrative of the conditions determining the different regimes that are possible 
in nonisentropic inflation models and it represents one of the main results in this paper. Figure [2] shows that there 
are no significant differences as regarding the dominating regions for each of the different polynomial forms for the 
inflaton potential, with little differences only between the cases p = 2, 4, 6. The figure also clearly shows that it is not 
always true that whenever T > H we are in a warm infiation regime of a thermal fluctuation dominated spectrum. 
The thermal bath can have a temperature much larger than the Hubble parameter and we still be in a quantum 
fluctuation regime when Q < 0.5. Likewise, the thermal bath can have in general a very low temperature compared 
to the Hubble parameter (or a temperature less than the one corresponding to the Hawking temperature) and the 
system can be in a warm inflation regime, where thermal fluctuations give the dominant contribution to the spectrum. 
These results come crucially from the fact that the quantum spectrum can depend on the temperature through a 
thermal enhancement factor, as shown in Eq. (j4.14p . which, as we have already discussed previously, comes from the 
consideration that in nonisentropic infiation models there is the presence of a radiation bath and that the inflaton 
fluctuations can be in thermal equilibrium with this produced radiation bath. 

B. Observational constraints to the power spectrum 

Having established in the previous subsection which are exactly the regions of parameter space where quantum or 
thermal fluctuations dominate, we turn now our attention to the question of the consistency of nonisentropic inflation 
in face of the recent observational data as coming from WMAP9yr. The observational parameters we test are the 
amplitude of the scalar curvature spectrum Aji{ko), the spectral index Ug and the tensor to scalar curvature spectrum 
ratio r, whose relevant expressions and data we have already given in section flV CI 

For the different forms of polynomial inflaton potentials considered, with powers p = 2 (quadratic), p — A (quartic) 
and p = 6 (sextic) that we have tested, we have produced plots for the two canonical values of number of e-folds 
before the end of inflation, A^* ^ 50 and N^, = 60. We have read (by numerically solving for) the spectrum of 
fluctuations leaving the horizon at these two values of N^. We have also produced sets of plots for the two cases 
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FIG. 3: The results for the V oc (f)^ potential: upper panels ((a) and (b)) have horizon exiting fixed at A^* = 50, lower panels ((c) 
and (d)), at A*'* = 60. In the panels in the first column ((a) and (c)), the dissipation term T is kept constant, while in the second 
column the dissipation is quadratically dependent on the inflaton field. In panels (a) and (b), Ai — 5.0 x 10^^^, A2 — 5.5 x 10^^^ 
and Acoid = 5.13 x 10"". In panels (c) and (d), Ai = 5.0 x 10"^^ A2 = 4.0 x 10"" and Acoid = 3.57 x 10"". The favorable 
regions by the observational data are the orange ones. 



of dissipation: a constant dissipation coefficient, T = Tg, and for a dissipation with a quadratic dependence on the 
inflaton field, T = C^cj)'^ /Alp, with the constant coefficients Tq and fixed by the numerical values taken for the 
ratio T/{3H) = Q. 

Instead of making the more traditional plots in a plane defined by the confidence levels (95% CL and 68% CL) for 
and r and specify the inflation models (for different parameters) as points in this plane, we have chosen, instead, to 
represent the allowed regions in the plane (Q, T/H). In the plane (Q, T/H), we plot for each of the cases of inflaton 
potential, TV* and dissipation coefficients, the permitted regions of parameter space consistent with the observational 
values of Ug and r. We have also draw lines of the inflaton potential parameter A in Eq. (|5.ip that are consistent with 
the amplitude of scalar perturbations A^{ko). Different color shades in the plots will represent the different regions 
of parameter space that are allowed or not. The obtained results for each case are shown in the Figs. [31 |3]and[5l in 
the case of including the thermal enhancement factor in the expression for r, Eq. (|4.29p . while in Figs. [51 [7] and [5] the 
thermal enhancement factor is absent. 

In all the Figs, shown in this subsection, we have used the following representation and notation. The gray and 
white areas in the plots will indicate the regions of parameters not favored by the observational data, while the 
different shades of orange will indicate those regions of parameters that are favored. The gray regions are the regions 
outside the lower and upper values (from the mean) for the spectral index, with value taken as Ug = 0.9636 ± 0.0084, 
obtained from the combined data coming from WMAP9yr+eCMB+BAO+i7o H- This value for Ug is within the 
68% CL region. The white areas are for r > 0.13 (corresponding to a region disfavored by the observational data at 
the 95% CL), for the cases of the quartic and sextic inflaton potentials, or for r > 0.08 (corresponding to a region 
disfavored by the observational data at the 68% CL), for the case of the quadratic inflaton potential. The different 
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FIG. 4: The results for the V oc cf>^ potential: upper panels ((a) and (b)) have horizon exiting fixed at A^* = 50 e-folds, lower panels 
((c) and (d)), at TV* = 60. In the panels in the first column ((a) and (c)), the dissipation term T is kept constant, while in the 
second column the dissipation is quadratically dependent on the inflaton field. In panels (a) and (b), Ai = 5.0 x 10"^^, A2 ~ 10^^^ 
e Acoid = 2.39 X 10~". In panels (c) and (d), Ai = 5.0 x 10"", A2 = 1.5 x 10"^^ and AcoW = 1.4 x 10~".The favorable regions 
by the observational data are the orange ones. 



shades of orange, from light to dark, separated by thin black lines, indicate the regions of parameters that have 
decreasing values for the tensor to scalar curvature spectrum ratio r. For the the quadratic inflaton potential (p = 2), 
we have 0.04 < r < 0.08 (light orange), 0.04 < r < 0.01 (medium orange) and r < 0.01 (dark orange), respectively. 
For the quartic {p = 4) and sextic (p = 6) inflaton potentials, we have 0.05 < r < 0.13 (light orange), 0.05 < r < 0.01 
(medium orange) and r < 0.01 (dark orange), respectively. Results for the amplitude of scalar perturbations are given 
by the long-dashed thick blue lines. These long-dashed thick blue lines are for three different values of the parameter 
A in Eq. (|5.ip that are consistent with the central value for the amplitude of scalar perturbations, = 2.41 x 10~^ 
(WMAP9yr-|-eCMB-|-BAO-|-i/o)- Acoid is the value required by the pure cold inflation case (Q = and T = 0), Ai and 
A2 are values below and above Acoid, respectively. All numerical values used for A are indicated in the figure captions 
for each case. Finally, the thick red line in each of the plots separate the regions dominated by quantum fluctuations 
(region below the line), or dominated by thermal fluctuations (region above the line), according to figured] 

In the figure[3]we have the results for the quadratic inflaton potential {p = 2). In the cold inflation regime, for both 
A^* = 50 and A^* = 60, this potential is already ruled out by observations (at 68% CL) Q in the pure cold inflation 
case (Q = 0, T = 0) and it is a viable potential in this case only at the 95% CL. However, for nonisentropic inflation, 
compatibility with observational data is recovered. For example, in the case of constant dissipation, compatible 
regions appear for low temperatures, T < H, and strong dissipation, Q ^ 1. In the case of inflaton fleld dependent 
dissipation, favored regions can run for all values of temperature, but for restrict values of dissipation, 0.1 < Q < 7.5 
(for A^* = 60), or 0.1 < Q < 2.5 (for A^* = 50). Regions of thermal fluctuation domination are more favored in the 
inflaton field dependent dissipation case only, while for the constant dissipation case quantum fluctuation tend to give 
the most dominant contribution only. 
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FIG. 5: The results for the V oc (fi^ potential: upper panels ((a) and (b)) have horizon exiting fixed at A^* — 50 e-folds, lower panels 
((c) and (d)), at A^, = 60. First column ((a) and (c)), the dissipation term T is kept constant, while in second column the dissipation 
is quadratically dependent on the field <f>. In panels (a) and (b), Ai = lO^^'', A2 = 10^^^ e Acoid = 3.59 x lO^^''. In panels (c) and 
(d), Ai = lO^^*", A2 — 10^^^ and Acoid = 1-78 x lO^^''. The favorable regions by the observational data are the orange ones. 



In the figure U we have the results for tlie quartic inflaton potential (p = 4). This potential is ruled out in the pure 
cold inflation case (Q = 0, T = 0), for both the = 50 and TV* = 60 cases, already at 95% CL when using the 
recent combined CMB data from from WMAP9yr+eCMB+BAO+i?o [ll- However, just like in the p = 2 case for 
nonisentropic inflation, this potential can also become again compatible with the observational data. All panels in 
figure 13] show compatible regions with quantum and thermal fluctuations domination. However, thermal fluctuations 
(warm inflation) dominated regions are more favorable in the field dependent dissipation cases, panels (b) and (d). 
The cases of constant dissipation, panels (a) and (c), show, however, much smaller compatible regions with the 
observational data and that it happens only in a thin strip region of Q values ranging from somewhere in the range 
1 < Q < 20, at least in the presence of the thermal enhancement factor in Eq. (|4.29p (see the results shown below, in 
the case when the enhancement factor in r is not included). 

In the figure [5] we have the results for the sextic infiaton potential (p = 6). This potential is completely ruled out in 
the pure cold inflation case {Q = 0, T = 0), for both the TV* = 50 and N^, = 60 cases, even when taking WMAP only 
results 0. This inflaton potential is the most constrained one by the observational data even in the nonisentropic 
inflation case. In particular, in the case of constant dissipation (panels (a) and (c) in flgure[5]), it is discarded by not 
having regions compatible with the upper bound on r, r < 0.13 at 95% CL (again, just like in the case of the quartic 
potential, this situation can change when discarding the enhancement factor in r, as we will discuss below). In the 
case of a dissipation coefficient that depends on the inflaton field (panels (b) and (d) in figure [5]), a compatible region 
with the observational data only appears for a small strip of values of Q, ranging in between 1 ^ Q ^ 10. Inside 
this range of Q, we find compatible regions for all values of temperature and also covering both quantum fluctuation 
dominated regions (that happens for T/H < 0.1) and thermal fiuctuation dominated regions (for T/H > 0.1). 

For the sake of comparison, in Figs. El E] and [5] we show the parameter space behavior when we do not take into 
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FIG. 6: Results for the <^ potential without the thermal enhancement factor in the tensor to scalar ratio r expression, otherwise, 
same notation and parameters used in figure [3] We can observe a much less restriction due to the tensor-to-scalar ratio (r) and in 
this case the range of Q is extended to small values. The favorable regions by the observational data are the orange ones. 



account the thermal enhancement term in Eq. (|4.29p . i.e., when considering a spectrum of gravitational waves that is 
decoupled from the thermal radiation bath generated during nonisentropic inflation. Note now that for all cases the 
region of parameters consistent with the observations increase, allowing as well larger regions for which r < 0.01 (dark 
orange shaded regions). These results show that the cases of constant dissipation tend to favor more the quantum 
fluctuation dominated regime, with very small or completely absent regions of thermal fluctuation dominated regime 
(like for a sextic inflaton potential, p = 6), while for a dissipation coefficient that depends on the inflaton field, we 
can have regions of parameters in both regimes. Note also that the larger region of favorable parameters appear when 
considering the quartic inflaton potential, with the sextic potential still being the most constrained case. 



C. Planck results 



As a final note, we notice that the recently released Planck results [s^ give for the tensor to scalar amplitude 
ratio the result r < 0.11 at 95% CL (when the high-£ CMB ACT-I-SPT data are added) and for the spectral index 
Hs = 0.9600 ± 0.0072, while when including the Planck lensing likelihood gives Ug = 0.9653 ± 0.0069 and r < 0.13, 
and by also adding BAO data, it gives Ug = 0.9643 ± 0.0059 and r < 0.12. These values are to be compared with the 
ones we have used here, Ug = 0.9636 ± 0.0084 and r < 0.13 (at 95% CL), obtained from the combined data coming 
from WMAP9yr-|-eCMB+BAO-(--ffo 01' which is quite consistent with those results from Planck. Therefore, all our 
results obtained previously remain fully valid when also accounting for the recent Planck data, indicating that the 
pressure on single field chaotic polynomial inflation models observed for cold inflation, does not hold in general for 
nonisentropic inflation models. 
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FIG. 7: Results for the (ji* potential without the thermal enhancement factor in the tensor to scalar ratio r expression, otherwise, 
same notation and parameters used in figure |4] The favorable regions by the observational data are the orange ones. 

VI. CONCLUSIONS 

With the increasing precision that are reaching the cosmological parameters obtained from the observational data, 
considerable constraints are been put on many inflationary models. With the recent release of the accumulated nine 
years of WMAP data, combined with other CMB data from ACT and SPT and other astrophysical data coming from 
BAO and precise measurements of Hq from large scale surveys, the simplest single field polynomial potential chaotic 
inflation models with quadratic and higher powers are ruled out at the 95% CL. Even simple monotonic inflation 
models appear disfavored at the 68% CL. These results, however, only apply to pure cold inflation models, in their 
most simple implementations and where the perturbation spectra is generated by quantum fluctuations only. In this 
work, we have analyzed a much larger class of inflation models, the noniscntropic inflation models. In these type of 
models the perturbation spectra is not sole determined by quantum fluctuations of the inflaton fleld, but because in 
this class of models a non negligible amount of radiation can be produced during inflation, thermal fluctuations can 
as well make a significant contribution to the power spectra. An example of a noniscntropic inflation model is warm 
inflation, where perturbation spectra come dominantly from thermal fiuctuations. Warm inflation models are known 
to have a much less constrained parameter space and these models can be made compatible with the observational 
data. However, no study has been made to date in order to explore parameter regimes ranging from pure cold inflation 
to warm inflation. We have made this analysis in this work. 

By appropriately making use of the stochastic inflation program, wc have successfully been able to describe both 
cold and warm inflation concomitantly. This has allowed us to obtain the combined contributions to the power 
spectrum coming from both quantum and thermal fluctuations. From these results, we were able to fully determine 
regions of parameter space when quantum fluctuations give the dominant contribution, or the opposite regime, when 
thermal fluctuations dominate the spectrum. 
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FIG. 8: Results for the (j) potential without the thermal enhancement factor in the tensor to scalar ratio r expression, otherwise, 
same notation and parameters used in figure [5] The favorable regions by the observational data are the orange ones. 



Since the observational data is strongly constraining and also putting considerable pressure on polynomial chaotic 
inflation type of models, which are also the type of models more easily prone to a microscopic quantum field theory 
description, as far model building for inflation is concerned, we have concentrated our analysis on these type of inflation 
models. We have, thus, investigated the viability of the noniscntropic polynomial chaotic inflation models in face of 
the most recent results for the cosmological parameters, including the amplitude of scalar curvature perturbations, 
the spectral index for the scalar perturbations and the ratio of tensor to scalar perturbations. We have investigated 
the cases of a quadratic, quartic and a sextic infiaton potential. In all cases we have found that the presence of a 
radiation bath and dissipation in the infiaton dynamics (which are the characteristics of all noniscntropic inflation 
models) can lead to compatible regions of parameter space with the observational data, re-establishing the viability 
of these infiaton potential models. 

Many characteristics of noniscntropic infiation models can be strongly dependent on the quantum field microphysics 
dynamics involved during infiation. These include, for example, the details of the model building, like how exactly the 
radiation bath is produced, the thermalization of the infiaton perturbations and the magnitude and dependence of the 
dissipation terms present during infiation. Even so, the approach we have adopted here is mostly model independent 
and our results can easily be extended to other type of infiation models, including, for example, hybrid infiation models, 
models motivated by string and superstring theories 0, [l^l, loop quantum gravity motivated infiation models 
braneworld type of models [52| - [53 |. infiation models with nonstandard potentials [55| . among many others. 
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Appendix A: Evaluation of the spectral index rig 



In order to determine the expression for ris, it is useful to define the following derivatives in terms of the slow-roll 
parameters: 
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where A = 4(1 + Q) + {Q — l)c. From the Friedmann equation for H, we obtain 
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Another useful expression is 
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For a given mode with wavenumber k, we can then write 



dk 
'dt 



k=aH 



Using also the equation for the entropy density in the slow-roll regime, 
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Ts = Q^'' , (A5) 
then, In r + In s = In Q + 2 In <j). Using also Eq. (jASp combined with the previous ones, we obtain the relation 
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Now we have to work out the dT /dk term. Assuming T = T(0, T), then 
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Using Eq. (jA6|) . we obtain: 
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Using c = T^tT/T, we obtain: 
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Using 
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(1 + 0) 



/3 and TH ~ 3H^Q, we obtain from the above equation, 
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If we want to evaluate Eq. (jA6p in terms of slow-roll parameters, c and dk/dt, Eq. (|A4p . we simply insert Eq. (|A10p 
in Eq. (|A6p . Note that for the cases of no temperature dependence on the dissipation coefficient, then we would have 
c = in the above expressions. Likewise, for a constant dissipation, wc would have /3 = 0. Analogously, we have: 
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